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^ Recently the identity method was proposed to calculate second moments of the multiplicity dis- 

tributions from event-by-event measurements in the presence of the effects of incomplete particle 

I 

identification. In this paper the method is extended for higher moments. The moments of smeared 
^ multiplicity distributions are calculated using single-particle identity variables. The problem of finding 
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the moments of the multiplicity distributions is reduced to solving of a system of linear equations. 
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I. INTRODUCTION 



The study of event-by-event fluctuations of chemical (particle-type) composition in high 
energy nucleus-nucleus collisions is a helpful tool to pin-down the properties of strongly in- 
teracting matter (see, e.g., review [T] and references therein). By measuring the fluctuations 
one may observe anomalies associated with transition between hadronic and partonic phases. 
Furthermore, the QCD critical point may be signalled by a characteristic pattern in the fluctu- 
ations [2HI]. Data on event-by-event chemical fluctuations from the CERN SPS [SHTj and BNL 
RHIC [8] were already published, and more systematic measurements are in progress. 

A serious experimental problem of event-by-event determination of the multiplicity of differ- 
ent hadron species is an incomplete particle identification caused by finite detector resolution. 
However, one can usually determine distributions Pj{x) of an identification observable x and 
obtain the average multiplicities (Nj) for different hadron species j. 

As an artefact of particle mis-identification measured fluctuations may be systematically 
reduced or enlarged. To overcome this problem a new experimental technique, called the 
identity method, was proposed in Ref. [H] for the analysis of events with two measured particle 
species. In Ref. this method was further developed to calculate the second moments of the 
multiplicity distributions of more than two particle species. In the present study we prove that 
the identity method can be generalized to determine third and higher moments of the multiplicity 
distributions in events consisting of an arbitrary number of different particle species. It was 
already emphasized that measurements of the third and higher moments for the event-by-event 
fluctuations is important for the search of critical point signals in nucleus-nucleus collisions 

mm- 

The paper is organized as follows. The used notations and definitions are introduced in 
Section [11} The identity method for the second moments of the multiplicity distribution is 



reviewed in Section |III[ In Section IV the method is generalized to third and higher moments 



of the multiplicity distribution. We consider in Section |V] the case in which the normalization 



of mass distributions are unknown, and finally in section VI we test the developed formalism 



in simulation. Section VII summarizes the paper. 



II. NOTATIONS AND DEFINITIONS 



We assume that particle identification (PID) is achieved by measuring a quantity Since 
any measurement is of finite resolution, we deal with continuous distributions of x denoted by 
Pj{x) {j = 1, . . . , k > 2, with k being the number of different particle types) and normalized as 

dx p,{x) = {N,) , (1) 

where (Nj) is the mean multiplicity of the particle type j. Note that the functions Pj{x) are 
found for the different particle species using the inclusive distribution of x for all particles in 
all events. The identity variables Wj{x) are defined as: 

»,(^) - ^ . (2) 

where 

k 

Pi^) = ^Pji^)- (3) 

i=i 

The complete identification (CI) of particles corresponds to distributions Pj{x), which do not 
overlap. This case leads to wj = for all particle species i j and wj = 1 for the j*^ species. 
When the distributions Pj{x) overlap, Wj{x) can take any real value between and 1. 
We further introduce the event quantities, Wj, for each particle type j: 

N{n) 

Wj ^ ' (4) 

i=l 

with j = 1, ■ ■ ■ , k , and define the event averages as 

{W,-...-W,) = J2w,-...-W,, (5) 

^^^^ n=l 

where N^v is the number of events and N{n) is the total multiplicity in the n*'' event. Each ex- 
perimental event is characterized by a set of measured values of PID variable x {xi, X2, . . . , xjv}, 
for which one can calculate the identity variables: {wj{xi),Wj{x2), ■ ■ ■ ,Wj{xN)}, with j = 

^ Here by x we refer to any experimentally measured particle property related to its identity, e.g. specific 
energy loss in a detector gas. 
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1, . . . ,k. Consequently, the quantities Wj, Wj, WpWq etc. are uniquely defined for each event, 
and their average values (see Eq. (|5])) can be found experimentally by straightforward averaging 
over all events. In the case of CI, one finds Wj = Nj, thus, Eq. ([s]) yields 

{W,-...-W,) = {N,.....N,) . (6) 

The ly-quantities are the event variables as they are calculated for each event. By averaging 
different combinations of these quantities over events, the moments of the ly-distributions 
{Wp . . . Wq) can be easily reconstructed directly from experimental data. However, the goal is 
to obtain the moments {Np . . . Ng) of the multiplicity distribution, which are unknown in the 
case of incomplete particle identification. The main idea here is to find the relation between the 
moments of the VT-quantities calculated from the measurements and the unknown moments of 
the multiplicity distribution. 

III. SECOND MOMENTS 

The relation between second moments of multiplicity distributions and corresponding second 
moments of Ws has been obtained in Ref. The quantities (Wj) and {WpWg) were found 
as specific linear combinations of all first and second moments (Nj) and (Nj), as well as the 
correlation terms {NpNg) . As the first moments are given by Eq. ([T]) the mathematical problem 
of finding of all the second moments (Nj) and (NpNq) is then reduced to solving a system of 
linear equations. 

IV. THIRD AND HIGHER MOMENTS 

In this section we derive formulae for all pure and mixed third moments of the unknown 
multiplicity distributions. Following the procedure developed for the second moments (see the 
previous section) we introduce the multiplicity {V{Ni, A^a,-)) ^^'^ ^^e PID probability density 
(Pj(x) = Pi{x)/ {Ni)) functions. Using these functions we get: 

Ni=oo N2=oo Nt^=oo „ „ 

Ni=0 N2=0 Nk=0 ^ 
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X J dxlP2{xl)... J dx%^P2{xU X ... X J dx^Pkixl)... X j dx%^Pk{x%^) 

[wp{x\) + ... + Wp{x\}^) + Wp{xl) + ... + Wp(a;^J + ^p(^i) + ••• + ^pI^at^)]^ 

k 

J]<iV/>^=^ + 3 J] {NlNi){w^;'-w^i)+2, J2 {N^Nl){w^;'-w^;) 



X 



1=1 l<i<KA; l<i<Kfc 



l<i<l<m<k i=l 

i<i<;<fc 

k 

-«v'-«w)+I]W(2«w' + '<-3<^-^ , (7) 



1=1 



X y" dxlP2{xl)... J dx%^P2{x%J X ... X j dx\Pk{x\)... j dx%^Pk{x%^) 

X [typ(xi) + ... + Wp{x\f^) + typ(xi) + ... + Wp{x\^) + Wp(xi) + ... + Wp{x%J\ 
X + ... + Wq{x\j^) + Wq{x\) + ... + + Wq{x\) + ... + Wq(a;^J] 

X + ... + Wr{x\f^) + + ... + w^(a;^J + ^^(a:^^) + ••• + Wr{x%y^ 

1=1 l<i<i<fe 

l<i<Z<fe 

l<i<i<m<fe 
fc 

i=l 

l<i<Z<fe 
k 

i=l 



{WlW,)= J2 E - E ^(^i,^2,...,iVfc) / dx\P,{x\)- dx\,P,{x\;) 
X J dxlP2{xj)... j dx%^P2{x%J X ... X j dx'lPkix^^)... 



X j dx''^^Pk{x%J [wp{x\) + ... + Wp(xjvJ + Wp{xl) + ... + Wp{x%J + Wpix'l) + ... + Wp{x%J] 

X [wgi^D + - + Wgi^lfJ + + - + Wgi^k) + ^^gi^l) + - + M^nJ] 

k 

= E(^'>^'-^+ E {N!N,){w^^.w^, + 2Wr)+ E {N^Nf){w^'.w- + 2WT) 



1=1 l<i<l<k l<i<l<k 



l<i<l<m<k i=l 



gp,i ' ^p,i 



l<i<l<m<k 1=1 

l<i<l<k ^ 
k 

+ E ^^'^ (2m^^ ■ ^ + - ^1,1 ■ ^ - 2«W ■ ^) • (9) 

i=l 

By exchanging in the last equation p by g one gets a similar expression for the {WpW"^"^. In 
Eqs. ([7]j9]) the following notations are used: 

= ■ Wj3 ■ , w'-p^^'' = ■ - ■ ■ , (10) 

1 f „ 1 

w7 



^ J w1{x)pj{x)dx , w'^q^i = J w^{x)wy{x)p,{x)dx , (11) 

w^i = J Wp{x)wg{x)wrix)pi{x)dx . (12) 

As it is evident from Eqs. ([7j|9| the third moments of the VT-quantities (M^p), {^p^g), 
{WgWp"^, and (VFpWqWr) are just linear combinations of the unknown third moments {N'p'^, 
(N^Ng), (N^Np), and (NpNgNr) of the multiphcity distributions. In Eqs. the first mo- 
ments (Nj) and second moments (Nj), (NpNg) are assumed to have already been determined, 



as it is discussed in Section |III[ The coefficients (10 12) entering into the right-hand-side of 
Eqs. ([^|9]) are also calculable, since they are expressed in terms of experimentally measured 
Pj(a;)-functions. The problem is then reduced to solving the system of linear equations for 
the third moments of the multiplicity distribution in terms of experimentally measured third 
moments of VT- quantities. This approach gives a unique answer, because by construction the 
number of equations is equal to the number of unknown third moments. 
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The method can be used to calculate any moments of the joint multiplicity distribution by 
deriving the corresponding set of linear equations. 



V. FIRST MOMENTS REEXAMINED 



As mentioned in sections [II] and III the first moments of multiplicity distributions can be 
directly obtained by integrating the corresponding x distributions, Pj{x). Here we consider 
the case when the shape of the Pj{x) functions are known, whereas their normalizations are 
arbitrary. This may be important for the analysis of hadron production in a narrow acceptance 
region. One then uses the data from the wider kinematic region to enlarge statistics and thus 
to improve the quality of the Pj{x) functions. We denote these arbitrary normalizations by Aj, 
i.e, 



dxp,{x) = {A,) . (13) 
It should be emphasized that all identities defined in Eq. ^ should be recalculated using the 



distributions defined in Eq. (13). 

Below we demonstrate that the identity method can also be applied to obtain the hadron 
multiplicities with these arbitrarily normalized mass distributions. The idea is to extract the 
unknown first moments of multiplicity distributions from first moments of ly-quantities. In- 
deed, one can write: 



W= E E - E 'PiNi,N2,...,N,) / dx\P^{xl)... / dx],P^{x\J 

Ni=0 N2=0 Nk=0 

X j dxlP2{xl)... j dx%^P2{x%J X ... X j dx\Pk{x\)... j dx)^^Pk{x%^) 

X [wp{x\) + ... + Wp(x]vJ + Wp{xl) + ... + + Wp{x\) + ... + Wp(x^j] 

k 

= 5](iV,)^, (14) 

where the mass probability distribution functions are now defined as Pi{x) = Pi{x)/ (Ai), and 

1 



Wp^i = I Wp{x)pi{x)dx . (15) 
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The equation above can be written in a matrix form: 



{W2) 



W2,l W2,2 ■■■ W2,k 



{N2) 



(16) 



\ i^k) J \Wk,l Wk,2 ••• Wk,k J \ {Nk) / 



The only unknowns in Eq. (16) are the first moments (Nj) which can easily be obtained by 
matrix inversion. 



Note that when {Ai) = {Ni) the right-hand-side of Eq. (14) yields: 



ppdx = (Np) . 



(17) 



i=l 1=1 

Therefore, in this special case the first moments of the multiplicity distributions are identical 
to the first moments of H^-quantities, (Wp) = {Np) , similar to the case of CI, cf. Eq. 



VI. TEST ON SIMULATED DATA 

In this section we apply formulae (7-9) to simulated data. As already mentioned in section [ll] 
we assume that particle identification is achieved by measuring a quantity x. In each event we 
generate three different particle specie^ called proton (p), kaon (K), and pion (tt). Further- 
more, the multiplicity distributions of these particles are simulated according to the Poisson 
distribution with different mean values. Consequently, the input information to our simulation 
is an inclusive (i.e., summed over all events and all tracks) distribution of the PID variable x 
for each particle (see shaded histograms in upper panel of Figure [T]) and mean multiplicities 
of particles. In our example these multiplicities are 8, 6 and 4 for protons, pions and kaons 
respectively. The simulation process comprises the following: (i) in each event we randomly 
generate from the Poisson distribution particle multiplicities using their mean values; (ii) we 
generate the quantity x for each particle j in an event from the inclusive distribution Pj{x) of 
X quantity for the corresponding particle. Using this simulated data we reconstruct all second 
and third moments of the multiplicity distributions. For each measurement x in an event we 

^ Note that the method does not have any hmitation on the number of particle species. 
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X [arb. units] 




FIG. 1. (Color online) Upper panel: The inclusive distribution Pj{x) of the particle identification 
variable x for protons, kaons and pions, used as an input to the simulation, are depicted with shaded 
histograms. The distributions of identities Wj{x) for protons, kaons and pions are presented with red 
solid, blue dashed and green dotted lines respectively. Lower panel: The distribution of event variables 
Wj with j standing for (a) protons, (b) kaons and (c) pions. 

calculate identity variables Wj{x) according to Eq. ([2]), where j stands for the particle type. 
Next, for each particle we determine the event variable Wj (see Eq. (ji])). Figure [l] (upper 
panel) shows the distributions of the identity variables Wj for different particle species. In the 
lower panel of Fig. [T] distributions of the Wj quantities are plotted. As seen from Fig. [T] the 
distribution exhibits evident structures. This is not surprising as by construction the latter 
is nothing else but the original multiplicity distribution smeared because of non-ideal particle 
identification, i.e., because of the overlapping Pj{x) distributions of different particles. These 
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structures are better seen for Wj^ because the x distribution of the pions is well separated from 
those of protons and kaons as illustrated by the shaded histograms in the upper panel of Fig. [T] 
Finally using Eqs. (7-9) we reconstruct all third moments of the particles from the corre- 
sponding moments of the W quantities. We would like to emphasize that all second moments 
needed for reconstruction of the third moments are obtained using the formulae of Ref. (TU] . The 
results are summarized in Table I, where the Poisson values are calculated using the moment 
generating function of the Poisson distribution, M{t) = e-<^>e<^>e*: 

< >= M"{t = 0) =< N >^ + < N >, (18) 

< A^^ >= M"'{t = 0) =< N >^ +3 < > + < N > . (19) 
In case of vanishing correlations between particles (as assumed in our simulation) one finds: 

< N^N2 >= (< iVi >2 + < Ni >) < N2 >, (20) 

< NiN2 >=< Ni >< N2 >, (21) 

< N1N2N3 >=< iVi >< N2 ><N-i> . (22) 

The statistical errors of the reconstructed moments of the multiplicity distributions result from 
the errors on the parameters of the fitted distributions Pj{x) of the identification observable 
X and the errors of the Wj quantitie^ Typically these two sources of errors are correlated. 
Moreover, one can experimentally study fluctuation signals built up from several moments of 
the multiplicity distributions. The standard error propagation is complicated and inconvenient. 
We therefore follow the Monte Carlo error calculation approach. We first randomly subdivide 
the simulated data into n subsamples and for each subsample reconstruct the moments M„ 
listed in Table I. In the second step we calculate the statistical error of each moment M as: 

a . 

(y<M> = (23) 



^ Contributions to the uncertainty of the reconstructed moments due to the assumptions on the functional 
representation of Pj{x) depend on the features of the experiment and their discussion is beyond the scope of 
this paper. 
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120.0275 ± 0.0320 


120.0174 
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192.0765 ± 0.0379 


192.0845 
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TABLE I. Results from the identity method for the 15-10^ simulated events with mean multiplicities 
corresponding to < p >= 8, < vr >= 6 and < K >= 4. All second moments (upper part) listed 
in the table are reconstructed according to Ref. |10j . all third moments (lower part) according to 



the formulae of Section IV The statistical errors were obtained from the Monte Carlo method using 
Eqs. (23-25). Moments directly calculated from the simulated multiplicity distributions are shown in 
the third column and the Poisson values obtained from Eqs. (18-22) are listed in the last column. 
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TABLE II. Results from the identity method for the 15-10^ simulated events with mean multiplicities 
corresponding to < p >= 8, < vr >= 6. In each event the number of kaons is taken to be equal to 
the number of protons in order to probe the method in presence of correlations. All second moments 
(upper part) listed in the table are reconstructed according to Ref. [10], all third moments (lower part) 



according to the formulae of Section IV The statistical errors were obtained from the Monte Carlo 
method using Eqs. (23-25). Moments directly calculated from the simulated multiplicity distributions 
are shown in the third column and the Poisson values are listed in the last column. Note that for 
protons and kaons Eqs. (20-22) are not satisfied anymore due to the introduced correlations. 
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where 

<M>=-y Mn, (24) 

and 
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FIG. 2. (Color online) Result of 100 simulations, each with 6-10^ events, with different mean multi- 
plicities of protons and pions while the multiplicity of kaons in each event is taken to be equal to the 
number of protons. In each box the ratio of the reconstructed third moment < > to the Poisson 



expectation (see Eq. (19)) is presented. 

The reconstructed moments, using the identity method, presented in Table [T] are in agreement 
with the generated moments within statistical uncertainties. Next, in order to demonstrate that 
the method works also in case of strong correlations, we introduced an artificial correlation 
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between protons and kaons. In each event we simulate the number of the protons and pions 
from the corresponding Poisson distributions as described above. However, the number of 
the kaons in each event is taken to be equal to the number of protons. The reconstructed 
moments from this simulation for 15-10^ events with the mean multiplicities corresponding 
to 8 and 6 for the protons and pions respectively are presented in Table [Tij Note that the 
mean multiplicity of kaons is the same as for protons as described above. As evident from 
the Table |TT] the reconstructed moments are again in reasonable agreement with the generated 
moments. However, for protons and kaons Eqs. (20-22) are not satisfied anymore due to the 
introduced correlations. 

In Fig. [2] we demonstrate the result of 100 simulations, with 6-10^ events in each, with 
different mean multiplicities of protons and pions. The number of the kaons in each event is 
taken to be equal to the number of protons, in order to induce a correlation between particles 
(see the text above). In each box the ratio of the reconstructed third moment < > of 



the proton multiplicity distribution to the Poisson expectation calculated using Eq. (19) is 
presented. One finds a high precision of the identity method independent of the values of 
average multiphcities. 



VII. SUMMARY 

In summary, we extended the identity method proposed in Refs. [9l [10] to third moments 
of unknown multiplicity distributions. We introduce the event quantities Wj depending on the 
single-particle identity variables Wj according to Eq. (|4|. Any moments of these quantities 
{Wp . . . Wg) can be measured experimentally. The problem of finding the second moments (Aj) 
and (NpNq) was solved in Ref. [lO]. The third moments of VT-quantities are found to be linear 
combinations of the unknown third moments (A'^), (A'^A^g) , (A^^Ap) and {NpNgNr). The 
first moments (A^) and second moments (Aj), (NpNg) entering into these relations are already 
known. We are confident that the procedure can be used for obtaining any higher moments of 
the multiplicity distributions by deriving the corresponding set of linear equations. We thus 
construct an iterative procedure in which all moments of the multiplicity distribution of the ra*'* 
order are calculated from measured n^^ order moments of the ly-quantities. By doing so, we get 
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a unique solution, as the number of linear equations in each step is always equal to the number 
of the unknown moments. We presented an explicit test of the method using simulated events 
including the case of strong correlations between particles. Note that the higher moments (e.g., 
third and fourth moments of pion and proton multiplicity distributions [HI [T2]) can reveal 
the effects of the QCD critical point with higher sensitivity in the event-by-event multiplicity 
fluctuations in nucleus-nucleus collisions 
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